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Abstract. Percolation and jamming phenomena are investigated for anisotropic sequential deposition of 
dimers (particles occupying two adjacent adsorption sites) on a square lattice. The influence of dimer 
alignment on the electrical conductivity was examined. The percolation threshold for deposition of dimers 
was lower than for deposition of monomers. Nevertheless, the problem belongs to the universality class of 
random percolation. The lowest percolation threshold (p c = 0.562) was observed for isotropic orientation 
of dimers. It was higher (p c = 0.586) in the case of dimers aligned strictly along one direction. The state of 
dimer orientation influenced the concentration dependence of electrical conductivity. The proposed model 
seems to be useful for description of the percolating properties of anisotropic conductors. 
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1 Introduction 

Physical properties of (partially) disordered systems are 
described in percolation approach [l|2|3|4|5|6|7j . Effect of 
concentration on physical properties is well understood. 
Nevertheless, alignment effect plays essential role for high 
aspect ratio objects, such as nanotubes and nanowires [8]. 

Percolating properties and behaviors of systems, com- 
posed of anisotropic nanoparticles, are extensively investi- 
gated during last years |9ll0lllll2ir3"IT3] . The problem of 
anisotropic percolation has been a subject of many inves- 
tigations. It is of interest and of value to inquire how an 
orientation of nanoparticles influences the main physical 
properties of the systems. Effect of the nanotube align- 
ment on conductivity is of specific interest. The anisotropy 
can be induced by various factors. For instance, alignment 
of the nanotubes may be induced by flow [8] or electric 
field pUE]. 

The electrical properties of a ceramic composition were 
considered in 12J. A profound influence of the conduct- 
ing phase structure on the electrical conductivity of the 
material was shown. The simulations indicated consider- 
able dependence of the percolation cluster properties on 
anisotropy of its components. 

The irreversible adsorption (deposition) of particles on 
solid surfaces is a subject of considerable practical im- 
portance. A well known example of an irreversible mono- 
layer deposition process is the random sequential adsorp- 
tion (RSA). RSA has attracted significant interest due to 
its importance in many physical, chemical, and biologi- 
cal processes. RSA is a natural model for irreversible and 



sequential deposition of macromolecules at solid-liquid in- 
terfaces. Some examples of the wide range of applicability 
of this model include adhesion of colloidal particles, as well 
as adsorption of proteins to solid surfaces, with relaxation 
times much longer than the deposit formation time. This 
process is well described in the literature and has been 
investigated extensively in the last decades. The topic has 
been well covered in |17|18j . 

Investigation of the irreversible adsorption of polyatomic 
species (A-mers) has received considerable attention in the 
last years. The results of the study of random sequential 
adsorption and percolation of polyatomic species on differ- 
ent substrates were presented in |13|19|20|21|22|23|24|25|26|27l28j . 
Numerical studies of random sequential adsorption (the 
RSA model) of rectangular particles on a flat substrate 
are performed and dependencies of the saturation concen- 
tration and the percolation threshold on the model pa- 
rameters are determined in |14j . 

The main goal of the present study is to investigate 
the influence of alignment on the percolation and jam- 
ming thresholds, as well as conductivity. In this paper, we 
provide accurate numerical data for the percolation and 
jamming thresholds of the (partially) ordered (aligned) 
dimers on a square lattice. We perform calculation of elec- 
trical conductivity as a function of dimers alignment. The 
proposed model can be useful for description of the per- 
colation behaviour of anisotropic conductor networks. 

The paper is organized as follows. In Section [2] the ba- 
sis of the model of deposition of dimers on a square lattice 
is presented, of The results, obtained using the finite size 
scaling theory, are also analysed and discussed in this Sec- 
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tion. The results, related to conductivity of the systems 
concerned, arc presented in Scction[3]pres. Finally, we dis- 
cuss dependence of the percolation threshold on the model 
parameters of interest in Section [4] 



2 Model and Simulations 

Different kinds of boundary conditions were used in perco- 
lation simulations: open, periodic or toroidal, and cylin- 
drical, i.e. open along one direction and periodic along 
another one. One uses usually the term crossing only to 
systems with open boundaries, and the terms spanning 
and wrapping for cylindrical boundary conditions [29j . In 
our study, the periodic boundary conditions are applied 
by gluing along the vertical and horizontal borders. We 
define a wrapping cluster as a cluster that winds around 
the system along the given direction, i.e. it provides a path 
with the length of 2tt [29] . 

Mersenne Twister random number generator [30] was 
utilized for filling in the lattice with dimer at given con- 
centration and orientation. It has a period of 2 19937 — 1. 

The anisotropy of dimer orientation is described by the 
orientation order parameter s: 



Ni + N- ' 



where Ni and 7V_ are the numbers of dimers, oriented in 
vertical and horizontal direction, respectively [31j . 

Let us consider a periodic square lattice of linear size 
L, on which dimers are deposited at random, but with 
given orientation. Two of nearest neighbour sites along 
given direction are randomly selected; if both sites are 
vacant, the dimer adsorbs on such sites. Otherwise, the 
attempt is rejected. In any case, the procedure is iterated 
until N ~ N\ + N- dimers are adsorbed and the desired 
concentration is reached. 

We perform our calculations of percolation thresholds 
using the Hoshen-Kopelman algorithm [52] . 

Fig. [1] demonstrates percolation in an anisotropic sys- 
tem of dimers (the dimers are oriented strictly along the 
vertical direction). 




The final state, generated by irreversible adsorption, is 
a disordered state (known as jamming state), in which no 
more objects can be deposited due to the absence of any 
free space of appropriate size and shape [IT]. If different 
orientations of the deposited objects are not equally prob- 
able, the definition of the jamming state is to be refined. 
Let us assume N\ > We define jamming for the fixed 
parameter s as a situation when there is no possibility of 
depositing any additional vertically oriented object. Nev- 
ertheless, there may be places for accepting horizontally 
oriented objects. 

Fig. [2] shows jamming for different s at a lattice of 
linear size L = 64. 





Fig. 2. Color online. Jamming at a lattice of linear size 64 at 
s = — 1 (left) and s — (right) (horizontally oriented dimers 
are shown in red, whereas the vertically oriented objects are 
shown in blue ) 



The percolation probability P vs. the site occupation 
probability p was obtained for the lattices of linear sizes 
L = 64, 128, 256, 512, 1024 and the number of runs 1000. 
The percolation probability P as a function of occupa- 
tion p for a particular value of s is shown in Fig. [3] The 
percolation threshold p c for a given lattice size L can be 
estimated from the condition P(p) = 0.5. 




Fig. 1. Color online. Percolation on a square lattice of linear 
size L — 64 at s — 1 (wrapping cluster is indicated in red) 



0,549 0,552 0,555 0,558 0,561 0,564 0,567 0,570 
P 

Fig. 3. Percolation probability P vs. occupation probability p 
for a lattice of linear size L = 1024 at s = 0.4, p c — 0.56141 ± 
0.00003 



The finite-size scaling analysis was carried out for get- 
ting the percolation threshold at L — » oo. The percolation 
threshold p c (L) was calculated for five different values of 
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the linear lattice size. The percolation threshold p c (oo) of 
an infinite lattice can be found by fitting the results for 
lattices of different sizes to the scaling relation (Fig. [4]): 



The average cluster size S and strength of the infinite 
network yield the scaling laws: 



|p e (i)-p c (oo)| cxL- 1 ^, 



(1) 



where the critical exponent v has the value 4/3 in three 
dimensions [BJ. 

In particular, our simulations gave the percolation thresh- 
old p c = 0.5863 for strictly oriented dimers (s = 1), and 
p c = 0.5617 for dimers oriented in two directions with 
equal probability (s = 0). The last result agrees with the 
known results |22I28| (see Table. [J). 



and 



S(p) oc \p-p c 



Pool?) °C {p~Pc) 



Table 1. Comparison of published and our results for dimers 
oriented in two directions with equal probability 





Pc 


Pjam 


Pc j Pjam 


Our result 


0.5617 


0.90668 


0.6195 


Open boundary condi- 
tions, L max = 2000 [22] 


0.562 


0.906 


0,62 ±0.01 


Periodic boundary con- 
ditions, L m ax = 112, 5 X 
10 4 runs [28] 


0.562 


0.907 


0.6196 




where 7 and /3 the universal critical exponents. In d 
2 E51: 



and 



7 = 2.3889, 
(3 = 0.1389. 



(4) 

(5) 
i = 

(6) 
(7) 



In our computations, the average cluster size and 
the strength of the infinite network demonstrate typ- 
ical behavior near the percolation transition. The critical 
exponents 7 = 2.02 ± 0.01, at p < p c , 7 = 2.37 ± 0.05, 
when p > p c , and (3 = 0.179 ± 0.002 for s = 0, extracted 
from the power laws ([4]) and are close to the known 
values © and 0. 

The fractal dimension of the incipient cluster [BJ, cal- 
culated as 







dr = d 

v 



, is close to the known values df — 1.896. For instance, 
d f = 1.8896 for s = . 

We found that parabolic law is a reasonable fit for p c 
vs. s (Fig. 03 . 

The jamming threshold for L — > 00 was found from 
scaling relation ([1]), where v — 1.0±0.1 is the critical ex- 
ponent [22j . In particular, for strictly oriented dimers (s = 
±1) (in fact, this case is quite equivalent to jamming in 
one dimension) Pjam 

= 0.8646, while exact value for jam- 
ming in one dimension is Pj am = 1 — e~ 2 ~ 0.86466 |17j . 
Fig. [5] demonstrates Pj a m(s). 



Fig. 4. Example of percolation threshold determination in the 
thermodynamical limit of (L — » 00) using the scaling rela- 
tion (© (s = 0, p c = 0.56173 ± 0.00003) 



Moreover, several quantities of interest [6] were calcu- 
lated. 

Average cluster size: 



s = y"v^ 2 /y^n^, 



(2) 



where rii is the average number of i clusters per lattice 
site. 

The strength of the infinite network P^, i.e. the proba- 
bility of an arbitrary site belonging to the infinite network: 



p n N °° 



(3) 





0,92- 




0,90- 


E 


0,88- 


CL~ 






0,86- 




0,84- 




0,59- 




0,58- 






CL 


0,57- 




0,56- 




0.55 -- 
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Fig. 5. Percolation threshold p c vs. orientation order param- 
eter s, p c (s) = 0.02442s 2 +0.56165, and jamming threshold vs. 
s, Pjam(s) = -0.0124s 2 - 0.0291|s| + 0.9065 



where N is the total number of sites. 
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3 Electrical conductivity 

The highly efficient algorithm, proposed by Frank and 
Lobb [33] was utilized for finding conductivity of a square 
lattice, filled with the dimers. 

The Frank and Lobb algorithm utilises the repeated 
application of a sequence of series, parallel and star-triangle 
(Y-A) transformations to the bonds of the lattice. The fi- 
nal result of this sequence of transformations is reduction 
of any finite portion of the lattice to a single bond that 
has the same conductance as the entire lattice. We used 
four equivalent resistors (conductors) with <Jf = 10 6 and 
<7, = 1 for occupied and empty sites [34], respectively, 
instead of each cell (see Fig. [6]). 



dir 


ler 






\» T 




















T 





Fig. 6. Equivalent schema of a square lattice filled with dimers, 
where 07 and <n are conductivity of the occupied side (07 = 
10 6 ) and empty one (at = 1), respectively 




Fig. 7. Conductivity a vs. occupation probability p for 
isotropic (s = 0)and strictly oriented dimers along the fixed 
direction (s — ±1). The lattice linear size is L — 256 



Examples of conductivity a as a function of occupation 
probability p at different values of parameter s are shown 
in Fig. [7J Behavior of the conductivity a(p) corresponds 
with direct estimations of the percolation threshold p c . 



There is a sharp transition of conductivity from <r_ to er + 
near the percolation threshold (see right inset in Fig. [7]). 
Moreover, a increases on passing from isotropic (s = 0) to 
strictly ordered dimers (s = ±1). 



10" 



10' 





+1 
-1 



10 10 10 

L 

Fig. 8. Application of scaling to calculation of the ratio of 
conductivity exponent to the exponent of the correlation length 
for isotropic (s = 0) and strictly oriented (s = ±1) dimers. 
Dashed line corresponds to the slope 0.973, which is typical 
for percolation in two dimensions 36 



It is known that conductivity a behaviour near the per- 
colation threshold (p — p c <C 1) obeys the scaling relations 



6j: 



cr_ oc (p c -p) Sc ,p < Pc, 

(7+ OC (p-Pcf C , P>Pc, 



(8) 
(9) 



where t c , s c are the conductivity exponents. 

A detailed study of the finite size effects is presented 
in order to discuss the universality class of the phase tran- 
sition which the system undergoes. The main aim of the 
paper is to determine the dependence 

For purposes of checking the universality class and cal- 
culating the values of t c ,s c , the relations JSJ) and ([9]) can 
be written in the following form 



er_ oc L Sc/v , p < p c , 



0+ 



oc L p > p c 



(10) 
(11) 



where L is the linear lattice size |36| . 

Fig. [5] shows results of the finite-size analyzes of er_ 
and er_|_. Ratios of the conductivity exponents to the cor- 
relation length exponent t c jv and s c /v are independent 
of s and close to 0.973 ± 0.05, which is a typical value for 
the percolation in d — 2 [36] . Hence, our model belongs to 
the class of random percolation in two dimensions. This 
statement agrees with the conclusion drawn from analysis 
of the exponents v , 7, (3. 

Near the percolation point, behaviour of the conduc- 
tivity depends essentially on s. At given p, the conductiv- 
ity can fall down near the percolation threshold (p — > p c ) 
if alignment of dimers (\s\ increases). Thus, alignment in- 
duces decrease of the conductivity at given conditions. 
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Fig. 9. Intrinsic conductivity k vs. parameter s. In the inset: 
(a/ffi — l)/p vs. p for s = 0, ±1 (symbols) and for conventional 
random site percolation (dash line). The data correspond to 
L — > oo. [55] , 

However, the situation changes drastically if occupa- 
tion probability is rather small p < 0.1. At given p, in- 
crease of s leads to decrease of a. The conductivity is 
maximum when all the dimers are aligned along the con- 
ductivity direction (see left inset in Fig. [7J. If p is small 
enough, the conductivity behaviour can be well fitted by 
virial expansion |37j 

a I a, = l + kp + mp 2 + (12) 

where k and m are the adjustable parameters. 

The value of intrinsic conductivity k can be extracted 
from (cr/ai — l)/p vs. p in the limit of p — > (see inset in 

Fig. ED. 

Intrinsic conductivity k as a function of orientation or- 
der parameter s is fitted well by linear function k = 1.88 + 
0.52s ( with corr elation p = 0.9976) (Fig. [9]). Notice that 
intrinsic conductivity calculated for conventional random 
site percolation (percolation of monomers) is 1.50 ± 0.01 
(dash line in inset in Fig. [9])). 



4 Conclusion 

Thus, new percolation problem, i.e. percolation of aligned 
dimers on a square lattice, was proposed and studied. The 
percolation threshold for deposition of dimers was lower 
than for deposition of monomer, p c — 0.5927. . ., never- 
theless, the problem belongs to the same universality class. 

The lowest percolation threshold p c = 0.562 corre- 
sponds to isotropic orientation of the dimers (s = 0). 
In the case of dimers aligned strictly along one direction 
(.s = +1,-1), the percolation point is p c — 0.586. Near 
the point of percolation transition (p — > p c ), the conduc- 
tivity essentially decreases if the absolute value of orienta- 
tion order parameter js| increases. Intrinsic conductivity 
k increases linearly with s (k — 1.88 + 0. 52s) and differs 
from the known value for a lattice filled with monomers 
k = 1.50 + 0.01. 



The proposed model can be applied to the phase tran- 
sitions of anisotropic objects on a lattice when their con- 
centration and orientation are varied. In particular, the 
model is useful for description of a phase transition from 
insulator to (semi)conductor upon aligned deposition of 
the prolate objects on a substrate. The natural extension 
of the model is substitution of dimers by fc-mers and inclu- 
sion of the bonds between the dimers into consideration. 



5 Acknowledgment 

This work has been supported by Russian Foundation for 
Basic Research (grant no. 09-02-90440) and Ministry of 
Education and Science of Ukraine (Project no F28. 2/058). 



References 

1. G. Grimmet, Percolation (Berlin: Springer- Verlag, 1999) 

2. H. Kesten, Percolation theory for mathematicians (Boston: 
Birkhauser, 1982) 

3. J. Feder, Fractals (New York: Plenum Press, 1988) 

4. J.M. Ziman, Models of disorder. The Theoretical Physics 
of Homogeneously Disordered Systems (Cambridge Univer- 
sity Press, 1979) 

5. B.I. Shklovskii, A.L. Efros, Electronic properties of doped 
semiconductors (Springer, Heidelberg, 1984) 

6. D. Stauffer, A. Aharony, Introduction to Percolation The- 
ory (Taylor & Francis, 1992) 

7. I. Sahimi, Application of Percolation Theory (London: 
Taylor &Francis, 1994) 

8. F. Du, J.E. Fischer, K.I. Winey, Phys. Rew. B 72, 121404 
(2005) 

9. D. L. Carroll, R. Czerw, S. Webster, Synthetic Metals 155, 
694 (2005) 

10. T.V. Sreekumar, T. Liu, S. Kumar, L.M. Ericson, 
R.H. Hauge, R.E. Smalley, Chem. Mater 15, 175 (2003) 

11. Y. Zhou, A. Gaur, S-H. Hur, C. Kocabas, M.A. Meitl, 
M. Shim, J. A. Rogers, Nano Letters 4, 1643 (2004) 

12. A.Yu. Dovzhenko, V.A. Bunin, Technical Physics 48, 123 
(2003) 

13. G. Kondrat, A. Pekalski, Phys. Rev. E 63, 051108 (2001) 

14. N.V. Vygornitskii, L.N. Lisetskii, N.I. Lebovka, Colloid 
Journal 69, 597 (2007) 

15. X. Liu, J.L. Spencer, A.B. Kaiser, W.M. Arnold, Current 
Applied Phys. 4, 125 (2004) 

16. C. Park, J. Wilkinson, S. Banda, Z. Ounaies, K.E. Wise, 
G. Sauti, P. T. Lillehei, J.S. Harrison, Journal of Polymer 
Science Part B: Polymer Physics 44, 1751 (2006). 

17. JW. Evans, Rev. Mod. Phys 65, 1281 (2003) 

18. Z. Adamczyk, Particles at Interfaces: Interactions, Depo- 
sition, Structure (Academic Press, 2006) 

19. H. Holloway, Phys. Rev. B 37, 874 (1988) 

20. JW. Evans, D.E. Sanders, Phys. Rev. B 39, 1587 (1989) 

21. Y. Leroyer, E. Pommiers, Phys. Rev. B 50, 2795 (1994) 

22. N. Vandewalle, S. Galam, M. Kramer, Eur. Phys. J. B 14, 
407 (2000) 

23. J. Cortes, E. Valencia, J. Colloid & Interface Sci. 252, 
256 (2002) 

24. V. Cornette, A. J. Ramirez-Pastor, F. Nieto, Physica A 
327, 71 (2003) 



6 



V.A. Cherkasova et al.: Percolation of the aligned dimers on a square lattice 



25. V. Cornette, A. J. Ramirez-Pastor, F. Nieto, Eur. Phys. 
J. B 36, 391 (2003) 

26. V. Cornette, A. J. Ramirez-Pastor, F. Nieto, Physics 
Letters A353,452 (2006) 

27. M. Quintana, I. Kornhauser, R. Lopez, A.J. Ramirez- 
Pastor, G. Zgrablich, Physica A 361 195 (2006). 

28. M. Dolz, F. Nieto, A.J. Ramirez-Pastor, Physica A 374 
239 (2007) 

29. G. Pruessner, M. Moloney, J. Phys. A 36, 11213 (2003) 

30. M. Matsumoto, ACM Trans, on Modeling and Computer 
Simulations 8, 3 (1998) 

31. D.A. Matoz-Fernandez, D.H. Linares, A.J. Ramirez- 
Pastor, EPL 82, 50007 (2008) 

32. J. Hoshen, R. Kopelman, Phys. Rev. B 14, 3438 (1976) 

33. D.J. Frank , C.J. Lobb, Phys. Rev. B 37, 302 (1988) 

34. Y. Yuge , K. Onizuka, J. Phys. C: Solid State Phys. 11, 
4095 (1978) 

35. A. Bunde, S. Havlin, Fractals and Disordered Systems, eds. 
A. Bunde and S. Havlin (Springer, 1996) 

36. C.J. Lobb , D.J. Frank, Phys. Rev. B 30, 4090 (1984) 

37. E.J. Garboczi , K.A. Snyder, J.F. Douglas, Phys. Rev. E 
52, 819 (1995) 



